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When a liquid droplet is located above a super-hydrophobic surface, it only barely touches the 
solid portion of the surface, and therefore slides very easily on it. More generally, super-hydrophobic 
q ' surfaces have been shown to lead to significant reduction of viscous friction in the laminar regime, 

so it is of interest to quantify their effective slipping properties as a function of their geometric 
characteristics. Most previous studies have considered flows bounded by arrays of either long grooves, 
or isolated solid pillars on an otherwise flat solid substrate, and for which therefore the surrounding 
air constitutes the continuous phase. Here we consider instead the case where the super-hydrophobic 
surface is made of isolated holes in an otherwise continuous no-slip surface, and specifically focus 
on the mesh-like geometry recently achieved experimentally. We present an analytical method to 
calculate the friction of such a surface in the case where the mesh is thin. The results for the 
effective slip length of the surface are computed, compared to simple estimates, and a practical fit 
is proposed displaying a logarithmic dependance on the area fraction of the solid surface. 
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Among the fascinating flow phenomena occurring on small scales ll|, |2j , super-hydrophobicity offers a unique bridge 
between microscopic features and macroscopic behavior 0, 0, 0, Q|. Super-hydrophobic surfaces are non-wetting 
surfaces which possess sufficiently-large geometrical roughness that a liquid droplet deposited on the surface would 
not fill the grooves of the surface roughness, but instead remain in a a fakir-like state where the droplet only touches 
the surface at the edge of the roughness (Fig. QJi). As a result, super-hydrophobic surfaces possess very high effective 



I. INTRODUCTION 
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. contact angles, and display remarkable macro-scale wetting properties 

One particularly interesting characteristic of super-hydrophobic surfaces is their low viscous friction. Since fluid 
... 

, in contact with the surface only barely touches it, but is instead mostly in contact with the surrounding air, small 
droplets can roll very easily, a phenomenon known as the lotus-leaf effect. In general, super-hydrophobic surfaces 
are expected to provide opportunities for significant drag reduction in the laminar regime, as has been confirmed by 
experiments cited below. The effective viscous friction of a solid surface is usually quantified by a so-called slip length, 
denoted here by A, which is the distance below the solid surface where the no-slip boundary would be satisfied if the 
flow field was linearly extrapolated, and the no-slip boundary condition corresponds to A = 00,S]- 

At low Reynolds number, the only characteristics affecting the friction of super-hydrophobic surfaces arise from 
their geometry, specifically (a) the distribution of liquid/solid and liquid/air contact at the edge of the roughness 
elements of the surface, and (b) the shape of the liquid/air free surface. For one-dimensional surfaces, i.e. surfaces 
which have one homogeneous direction (Fig. [TJd), significant drag reduction can be obtained when the homogeneous 
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mogene ous 

3H,I3 



direction is parallel or perpendicular to the direction of flow, as demonstrated experimentally 

and theoretically 0, H 0, 0, 0, 0, H- 

For two-dimensional surfaces, an additional free parameter is the topology of the air/solid partition on the planar 
surface, and two general types can be distinguished. In the first type, the air is the continuous phase, and the liquid 
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FIG. 1: (color online) Typical topology of super-hydrophobic surfaces: (a) Side view of a random super-hydrophobic surface, 
illustrating the contact of the liquid (gray) at the edge of the surface roughness (blue); (b) Top view of a one-dimensional 
surface; (c) Top view of regular two-dimensional surfaces where the air (white) is the continuous phase; (d) Top view of regular 
two-dimensional surfaces where the solid (blue) is the continuous phase. The case illustrated on Fig. QJl (top) is the focus of 
the current paper. 



is in contact with the solid only at isolated, unconnected, locations (Fig. \Vp). This is the most commonly-studied 
type of super- hy drop hobic surf ace, and arises for surface roughness in the shape of bumps or posts on an otherwise 



flat solid surface [lfj, | 
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The second type of surface, less studied, is one where the solid is the 



continuous phase, and the liquid/air contact occurs on isolated domains (Fig. [T]i). These surfaces can be obtained 
by making holes in an otherwise flat material [2^ , and have been used to study the influence of the geometry of the 
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liquid/air free interface on the viscous friction of the surface 

A super-hydrophobic surface with a continuous solid phase can also be obtained by using an intertwined solid 
mesh. Recently this method has been exploited experimentally, using coated steel 0, [31] , and coated copper [3^] , as 
reproduced in Fig. [2K and b. The resulting surface is also the relevant geometrical configuration for flow past textile 
and fabric material. In this paper, we present an analytical method to calculate the friction of such a mesh-like (or 
grid-like) surface in the case where the mesh is thin in the plane (mesh aspect ratio e <C 1). The analysis, based on a 
distribution of flow singularities, leads to an infinite system of linear equations for the approximate Fourier coefficients 
of the flow around the mesh, and gives results with relative error of order eln(l/e). The results for the effective slip 
length of the surface are computed, compared to simple estimates, and used to propose a practical fit displaying the 
expected logarithmic dependance on the solid area fraction 
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II. SHEAR FLOW PAST A RECTANGULAR GRID 



A. Calculation background 



Our paper follows classical work on the broadside motion of rigid bodies at low Reynolds number in the asymptotic 
limit for narrow cross-sections 
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34j. Here, we build on the method of Lcppington and Levine who studied axisym- 



metric potential problems involving an annular disk by using a distribution of singularities, and obtained an efficient 
method in the limit where the radii difference of the disk approaches zero 35|. Roger and Hussey studied the flat 
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FIG. 2: (color online) Mesh-like super-hydrophobic surface, (a): Coated copper mesh at the micron scale [32|; (b): The mesh 
from Figs. [2^i is super-hydrophobic (water droplet, volume about 4 fjl) [3^]; (c): For our calculation we consider a planar square 
mesh of periodicity b and width eb; (d) : A shear flow is set up in one of the principal directions of the mesh, with shear rate 7; 
the far-field flow is given by v = j(z + X)e x , where A is the effective slip length of the mesh. Figs. [2^, and [2b reprinted from S. 
Wang, Y. Song and L. Jiang (2007) "Microscale and nanoscale hierarchical structured mesh films with super-hydrophobic and 
superoleophilic properties induced by long-chain fatty acids", Nanotechnology, 18, 015103, courtesy of Shutao Wang and with 
permission from IOP Publishing. 



annular ring problem both experimentally and theoretically using a beads-on-a-shcll model to represent a distribution 



of point forces [36l | . Subsequently, Stewartson showed that viscous fluid exhibits a marked reluctance to flow through 
a thin torus and the flux function has an essential singularity in the limit in which the hollow boundary approaches 
a circle and disappears [37[. This 'blockage' feature was further demonstrated by Davis and James in considering 
flow through a array of narrow annular disks placed in planes normal to the flow in order to model a fibrous medium 
j^ . Their analysis was made tractable by retaining only the inverse square root term in the force density function. 
A similar asymptotic estimate was used by Davis in modeling the broadside oscillations of a thin grid of the type 
discussed below [3^]. Since edgewise motions induce the same edge singularity in the force density, we follow in this 
paper a similar analytical treatment. 



B. Model problem 



We consider a thin stationary planar square mesh in the [x, y) plane subject to a shear flow, with shear rate 7, 
along the x direction, which is one of the principal directions of the mesh, as illustrated in Figs. and d. The mesh 
has periodicity b and a small width eb. We denote by e x and e y the unit vectors parallel to the mesh axes, and e 2 is 
the third Cartesian vector. We further assume that the shape of the interface between the liquid and the air located 
underneath the mesh is planar, with no protrusion of the mesh inside the fluid. Since the mesh is a super-hydrophobic 
surface and presents to the flow a combination of no-slip domains (the mesh) and perfectly-slipping domains (the 
free surfaces), the effective flow past the surface can be described as a slipping flow, and is given in the far field by 
v = A f(z + X)e x , where A is the effective surface slip length (A > 0). The goal of the calculation below is to present an 
accurate analytic calculation for A in the limit where e«l. 

In order to perform the calculation, we move in the frame translating at the steady velocity Ue Xl with U = 7A. In 
that case, the problem is equivalent to the uniform translation of an thin mesh in its own plane at velocity — Ue Xl 
and the goal is to derive the value of the shear rate 7 in the far-field, \z\ 3> b. The slip length will then be given by 
A = U/'j. Edge effects are ignored by assuming an infinite mesh as then periodic point force singularities can be used 
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FIG. 3: (color online) A typical grid element (see text for details). 



to describe the fluid flow field. The rectangular grid naturally introduces Fourier series with respect to Cartesian 
coordinates, and the analysis establishes a set of integral equations with the same logarithmic dependence on the 
geometrical parameters of the grid as in the flows described in the references discussed in ^11 Al 

The Reynolds number of the viscous incompressible flow is assumed to be sufficiently small for the velocity field v 
to satisfy the creeping flow (Stokes) equations 33, 3| 



^V 2 v = Vp, V • v = 0, (1) 

where (i is the coefficient of viscosity and p is the dynamic pressure. The prescribed grid velocity is given by 

v = -Ue x at z = 0, (x, y) e G. (2) 

In a typical grid element < x, y < b, G is complementary to the square hole bounded by the lines x or y = \eb or 
(1 — \t)b (shown in black in Fig. [3^i). 



C. Solution using a superposition of singularities 



Following the calculation for broadside oscillations of the grid in Ref. [39(, the fluid motion can be represented as 
due to a distribution of tangcntially directed Stokcslets over the flat mesh G and the density functions must be both 
periodic in two dimensions and symmetric with respect to the sides of each square. The field due to a two-dimensional 
square array, period b, of point forces of strength AirfiUb directed parallel to the mesh's motion, is governed by 



V • v A = 
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and 
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//V 2 v^ — — 4:7TfiUbe x S(z) S(x — nib)5(y — 77,26) 
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= — 7— e x 6(z) 2^ exp 



mi— — OO 7712 — — 00 

With m = m\G x + m2e y , the m = term in Eqs. ([3]) and (H]) yields 

Att/iU 



2m. 

— [mix + m 2 y) 



e x 5(z), V-v o = 0, 



whose solution, 



2ttU\z\ 

v = 7 — e x , pq = 0, 



(4) 



(5) 



(6) 



exhibits the anticipated shear at infinity. Note that the suppression of the immaterial arbitrary multiple of Ue x in v 
ensures uniqueness below. The solution of 



2uU f°° ■> 

AtV 2 v - Vp = —r — e x / e lkz dk ^ ' exp 

" J — OO m 



2m 



-m • r 



V-v = 0, 



with a prime denoting that the m = term is omitted, is readily found by Fourier transform techniques (see Ref. 40j ) . 
Thus the flow governed by Eqs. © and (01 is compactly expressed as 



v A = U 



2n\z\ 



Si)e x + V 



.dS2 

dx 



PA = HU—, 
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Only the velocities at the mesh and at infinity are needed for the subsequent analysis. The solution Eq. ([8]) shows 
that 



2nU\z\ 



-e x = v as \z\ — > oo, 



d 2 S 2 
dx 2 



since all terms in Eq. © exhibit exponential decay, and 

[va] z =o = Ue x (-Si 
where, after substitution of Eqs. ([9]) and (fT0|) , we have 

C m = 1 



-Ue x V'- — -exp 
^ m 



2?ri 



C m , 



2 m : 
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In particular, we see that 



C(mi,0) — 2' ^ n 



as mi — > cx>, 



C, 



(0,m 2 ) 



1, Cm — > 1 as 7712 — > OO. 



(14a) 
(14b) 



The shaded region in Fig. [3t> suffices for the distribution of periodic point forces and the enforcement of the 
prescribed mesh velocity (Eq.EJ). The hexagon in Fig. [3b lying along the x-axis is given by \y\ < ieb, \y\ < x < b — \y\, 
or — 1 < w, s < 1 in terms of new variables (w, s) defined by 



-6 — x 1 1 

y = 2 ebs > x = ~ ^ ~ es ^ 



(15) 



Here w may be identified as the tangent of an angle subtended at the center of the square. The hexagon in Fig. [SJd 
lying along the y-axis is given similarly by interchanging x and y in Eq. (JTSj). The flow generated by the translating 
mesh can be written as 



1 



f x (w, s)v A 



'-l J-i 

+ fy{w,s)VA 



x — 2^(1 — w + es w), y — 2 £ ^ s ' z 



x — 2 6 ^ s ' y ~ — w + esw), z 



(16) 



where f x and f y are dimensionless force densities. The prescribed mesh velocity is then obtained by enforcing Eq. © 
at each point of the two hexagons. Thus Eq. (fT6|) gives 
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for all — 1 < W, S < 1. We further note that the summation in Eq. (fl~2"|) takes the simpler form 
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from which we note that large contributions to Eqs. (JTTJ) arise from summations of type 
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Substitution of Eq. IT2|) into Eqs. (fT7|) and use of Eqs. (fT4]) now yields, when terms that tend to zero as e — > are 
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for all — 1 < W 1 S < 1. All the w-intcgrals in these integral equations can be expressed in terms of Fourier coefficients 
defined as 



f 



lfxn(s), fy„(s)} = - J [f x (w, s), f y (w, s)] cosriTriuw (n>0), 



(21) 



because the symmetric forcing ensures that the corresponding sine coefficients are all zero, and where we use brackets 
to define simultaneously two sets of Fourier coefficients. We then identify Eqs. ([20]) as a pair of Fourier cosine scries 
in W, for each S in [—1, 1], and, by considering coefficients of C0S7ttoiW and cos Trm2W in the respective series, we 
obtain 
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for all —1 < 5 < 1. These two related sets of equations have the form 
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whose solution is 



with the particular property 



[/»»(«), fyn(s)} = ^ n jf 2 f 1/2 ("!<«< 1), 



[fxn{s), fyn(s)]s = [A n , B n }. 



(24) 
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Thus, as demonstrated more rigorously by Lcppington and Levine [351 ] and exploited frequently elsewhere, the loga- 
rithmic kernel obtained as an asymptotic estimate yields the inverse square root function as the associated asymptotic 
solution. The substitution of Eqs. (|23|) and (|25|) into Eqs. (|22| finally yields an infinite system of linear equations 
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D. Determination of the slip length 



The flow at infinity is determined by substitution of Eq. (fTT]) into Eq. (fl6|) , which gives 



2irU\z\ 1 

v ~ ; e x - 

b 4 



[f x (w, s) + f y (w, s)]ws as |z| 



(27) 



whose simple form is due to the exponential decay of all Fourier modes except m = 0. Then Eqs. ([2Tj) and (|25|) show 
that the shear rate at infinity is given by 



7 : 



7rC/(A + B ) 



and therefore the slip length, A, is found to be 



7 7r(,4 + B ) ' 



(28) 



(29) 



which is proportional to the only other length scale in the problem, namely b. Notably, only the zeroth-order coefficients 
of Eqs. (f2"6"| contribute to the slip length. The system given by Eqs. (f2"6"|) . independent of b, is solved in truncated 
form for various values of e in EjTTTJ 
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FIG. 4: (color online) (a): Dimensionless effective slip length for the mesh-like super-hydrophobic surface, X/b, as a function of 
the aspect ratio of the mesh, e. (b): Ratio between the slip length given by solving Eqs. (|26|l and the simple estimates provided 
by Eq. ((33]) (A/Ai, blue, solid line) and Eq. (f36|) (A/A2, black, dash-dotted line), as a function of the aspect ratio. 

E. Error estimate 

An error estimate can be gleaned by retaining all terms in proceeding from Eqs. (|17|) to Eqs. (f2"2"| . Although the 
former suggests that the symmetry of each grid element implies a mathematically even dependence on e, convergence 
considerations prevent the error bound from being 0(e 2 ). For example, the Fourier cosine series for (ir — x) 2 in [0, 2tt\ 
shows that 



5—5 = o mm ( s ) - emax(s, 5) = 0(e . (30) 

m— 1 

Hence the (relative) error estimate is 0(eln(l/e)), as in various sample calculations given by Leppington and Levine 



351 ]. Mathematically, its presence is due to an elliptic integral generated by the fundamental singularity. 



III. NUMERICAL CALCULATION AND ANALYTICAL ESTIMATE OF THE MESH SLIP LENGTH 

A. Numerical results 

We now solve numerically the infinite series given by Eqs. (|26p by truncating it at a finite value m\ = m-2 = N. 
The numerical results are displayed in Fig. [3^, where we plot the dimensionless slip length, A/6, as a function of the 
aspect ratio of the mesh, e, for a truncation size of N = 10 2 . The exact size of the truncation has little influence on 
the final computed results; for example, when e = 1CP 3 , the relative change in the computed slip length is less than 
0.01 % between N = 10 2 and N = 10 3 . 

Our results show that the surface does show a reduction in friction (i.e. A > 0), and the effective slip length, always 
on the order of (but typically smaller than) the mesh size, increases when the solid area fraction decreases. Since the 
data on Fig. 0^ are plotted on a semilog scale, we see that we recover the expected logarithmic relationship between 
the effective slip length of the mesh and its aspect ratio A least squares fit to the data shown on 

Fig. leads to the empirical formula 

^■ = Alne + n, (31) 


where the values A = —0.107 and II = —0.069, gives result with maximum relative error of 0.4%. Alternatively, if we 
define the solid area fraction of the surface <fi s , it is easy to see that <f> s = e(2 — e). A least-square fit of the type 

^=A 2 ln0 s +n 2 , (32) 
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with the parameters A2 = —0.107 and II2 — 0.003, leads to a maximum relative error of 0.6%. As a matter of 
comparison, the fit proposed in Ref. [261 ]. based on numerical simulations for <p s > 0.05, leads to a relative error of 
5% with the small-^ results obtained here; the asymptotic calculations presented in this paper are therefore able to 
quantitatively agree with the results of Ref. [26j which are valid for much larger mesh sizes. 



B. Comparison with simple estimates 

An estimate of the mesh slip length can be obtained analytically by performing the truncation by hand. If only 
Aq, Bq — that is the mean values of the density functions with respect tow — are retained in Eqs. (|26p . the surviving 
two equations give the lowest-order estimate 

f = n(A + B a ) = 3^ ln (^)' (33) 

We define here L = ln(2/7re) for subsequent convenience. Note that the slope of the relationship given by Eq. (|3"3"|) is 
1/3tt « 0.106, which is very close to the fitted value |A| = 0.107 in Eq. (pH"j) . 
Further, on writing Eqs. (|26a|) and (|26c[) as 

/ 9 \ 00 f_1 )m 

3 = (A + B ) In — + V [ -^(A m + B m ), (34a) 

V ' m=l 

/ 2 \ 00 f_i yn 

-1 = (Ao -Bo) In — - V ^-{A m - B m ), (34b) 



m—1 



and observing that the matrix elements in Eqs. f)26j) have L only on the diagonal, it may be deduced that 

y-^ = o(i" 1 ). (35) 

As an illustration of this result, the truncation of Eqs. (|26[) at mi = rri2 = 1 shows that the leading-order correction 
to Eq. (|33|) leads to the improved estimate 

A2 - L - - (36) 



b 3vr 6 ^ 



On Fig. HJd, we show the ratio between the actual slip length (as calculated in Fig.[3Ji), and the simple estimates from 
Eq. (f3"3"| (blue, solid line) and Eq. (f3"6"| (black, dash-dotted line). We see that these simple formulae over-estimate 
the slip length (by up to 10%), but asymptotically converge to the correct result when e — > 0. As expected, the 
second-order solution, Eq. (|36[) . is quantitatively better than the first-order one, Eq. (l33l) . 



IV. CONCLUSION 

In this work we have considered super-hydrophobic surface in the shape of a square mesh, and have presented an 
analytical method to calculate their effective slip length. Our analysis, which is valid for asymptotically small aspect 



ratio of the mesh, agrees also quantitatively with numerical calculations valid for larger mesh width [26j . The results 
we obtain show that, for thin meshes, such surfaces can provide slip lengths on the order of the mesh size. However, for 
realistic mesh aspect ratio, say 10% or larger, the slip length is only a fraction of the typical mesh length (see Fig. [4|. 
The increase of the slip length with the solid-to-air ratio of the overall surface is logarithmic, a universal feature of 
super-hydrophobic coatings with long and thin no-slip domains 

0, EE 

l22j. Beyond such physical scaling, our 



results provide a first analytical prediction for the resistance to shear flow of the recently devised super-hydrophobic 
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surfaces made of metal wires [5|, |3l| , and more generally to geometries with grid- like features such as fabric and textiles 
composed of interwoven threads. 

Although we have considered here an idealized model system, the analytical approach developed in the paper has 
the advantage that it provides the effective slip length of the surface with relative error of order elog(l/e), which is 
significantly better than the more traditional but only logarithmically correct slender-body approach. 
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